VOL. 20, NM" 1 JUMHO 2004 159M 0717-2915 


A] => 






Jour 


A IR 
A A A A 
BOGIETT FOR REMOTE RENBINO AMD SPACE MPORMATION SYETEMS 









a e -M 
¿E ha A 








mí DÁ ll EL MN 


CAPITULOS CONSTITUIDOS DE SELPER 





E li 





CAPITU LOS 
ESPTCIALIS 








CAPITULOS EN 
FORNACIÓN 





Ml BM 


SOCEDAD LATINOAMERICANA DE PERCEPCIÓN REMOTA 
Y SISTEMAS DE INFORMACIÓN ESPACIAL 


SOCEDADE LATINO AMEFICANA FM SENSORIAMENTO REMOTO 
FE SISTEMAS OF IMPORMACADO USPACIAS 


LATMAMERFICASN SOCIETY FOR REMOTE SIONIMG 
ANO SPACE INFORMATION SYSTEMS 


AGRADECIMIENTOS - AKNOWLEDGEMENTS 


as kl 







EN 2 


| ES y ESA - Asuntos rminrnacionales 
8-10 rue Mario Nihis 





ESA +» ESHUIM 
Vis Galileo Galde! 
Casañas Postañe LA 


AY 
"EN 
00044 Fresco (Momajl 
E ITALIA 


Tel: (19 05) 944180 940 


Fax: (15 041 14110 220 
ll es al Att: Se Maurino Foa ¡Armáérica Latina! 


y 





157318 PAÑOS Codex 15, FRANCIA 
Tel: (33.11 5300.T7T31, Fax: (19.1) 1909 7027 
Ar Sra. Mario E iabeth De Vol [América Latina) 





REVISTA - JOURNAL SELPER VOL. 20, N* 1, 2004 








REVISTA. JOURNAL 


VOL. 20 WN 1 ¿JO 7004 5 0717.3014 











q” PORTADA 


Macizo de San Rafael 


Subescena ASTER RGB:321 de 15 m de 
rosolwción espacial. Se observa un sector 
del Macizo de San Rafael ubicado a unos 
40 km al veste de la ciudad de San Rafael, 
Provwncia de Mendoza, caracterizado por 
depósitos del Proterozorco Superior y 
Paleozoico. Ál suroeste se encuentra la 
estructura concéninca del pórtico de cobre 
y molibdeno conocido como El Infernllo. Al 
norte se puede apreciar el dique Los 
Reyunos (34* 36' 8 18 Lattud Sur y 68" 35" 
31.27" Longitud Cesto) que ombaisa las 
agues del Rio Diamante 
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A parir de las deasiones adoptadas en el marco del X Simpo- 
sio Latinoarrercano de Perceación Remola y Sistemas de Intor- 
mación Espacial, organizado por nuestra Saciedad; llevado a 
cabo en Cochabamba, Bolesa, en noviembre de 2302, al Comite 
Ecitcris! ha tomado el compromiso de regularizar las publicacio- 
nes de la Revista Técnica Clentífica. 

En este sentido es necesarío dasacar el invalorable apoyo 
oblenido de la Agencia Espacial Eurcoea ¡ESA), quien ha toma- 
do baje su responsablilcad ta impresión de los volúmenes 16, 
17, 16 y 19 correspondientes a los años 2000, 2301 y 2002 
respecilvamente los cuales se encuentran disponibles para los 
miembros de nuestra Sociedad. 


En esta oportunidad hacemos begar la publicación de volumen 
20 numero 1, dorde se incluyen babajos que han sido erviados 
por sepecalistas en el 3rea de Teledetección y Sistemas dea Infor- 
mación Geográfica que desarrollan sue actividades en diferentes 
institaciones de América Latina, las cuales han sido remitidos a 
expertos evaluadores para su ravisión y aprobación, Nos loca 
agradecer muv especialmente la labor da estos expertos que 
posibilitan un mejorarientc en la calidad de os articulos que en 
ella se presentan. 


Nos toca agradecer muy especialmante la lasor llevada a cabo 
por el grupo de evaluadores que ha participado en la revisión de 
105 MISMOS. 
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Los articulos recibidos serán enviados s tres (34 
expertos en la temática para su revisión. Los tra- 
bajos aprobados serán publicados en estricto 
orden, de acuerdo a las fechas de llegada de las 
contribuciones. 

Los idiomas oficiales SELPER son: Español, 
Portugués e Inglés. 

Los trabajos deberán estructurarse contem- 
plando las siguientes secciones: 


a) Título del trabajo. Nombre de los autores y 
direcciones completas 

b) Resumen (no más de 150 palabras) indican- 
do al final las palabras claves. Deberá incluirse en 
Español o Portugués, además de Inglés 

c) Introducción 

d) Objetivos 

e) Metodología empleada y materiales 


fi Resultados obtenidos 

g) Conclusiones 

h) Bibliografía: sólo se incluirá la citada en el tex- 
to. Se indicarán los autores, por orden alfabético, 
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incluyendo volumen y páginas, cuando corre- 


sponda. 
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lo y estar numeradas en forma consecutiva. 

Se deberá enviar una copia del trabajo en for- 
mato Word y una copia papel. La extensión total 
del trabajo no deberá superar las 12 páginas, 
(DIN-A4). 
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Comparativo de fusión de datos para la 
elaboración de cartografía predial urbana 


Judith Ley Garcia 


Laboratorio de Geomática, Instituto de Investigaciones Sociales, Universidad Autónoma de Baja 
California, Edificio de Investigación y Posgrado. Blvd. Benito Juárez s/n. C.P, 21280. Mexicali, B.C., 
México. Tel:(686)-566-29-85 Ext: 130. Fax Ext : 135, E-mail: ¡leyívuabc.mx 


RESUMEN 


Las técnicas de fusión de datos, brindan la posi- 
bilidad de integrar imágenes con diferentes carac- 
terísticas espaciales y espectrales, para obtener 


explorar tres técnicas básicas de fusión de datos 
(multiplicativa, sustitución del componente princi- 
pal, sustitución del canal de luminosidad) para inte- 
grar las imágenes pancromática y multiespectral 
del sensor Quickbird; y comparar el potencial de 
las imágenes resultantes para su utilización en la 
elaboración de cartografía urbana a escala predio 
de la ciudad de Mexicali, Baja California, México. 


Palabras clave: Fusión de Datos, QUICKBIRD, 


Percepción Remota, Cartografía Urbana, Monito- 
reo Urbano, Nitidez, Resolución Espacial. 


The data fusion techniques, give the possibility of 
integrating images with different space and spectral 
characteristics, to obtain secondary images of 
better visual quality than the original ones. The pres- 
ent work is centered in exploring three techniques 
basic of data fusion (multiplicativa, 
the principal t, substitution of the lumi- 
nosity channel) to a po the images 

tica and Il of the Quickbird 
sensor; and to compare the potential of the resulting 
images for ¡ts use in the elaboration of urban cartog- 
raphy of the city of Mexicali, Baja California, Mexico. 


“Keywords: Data Fusion, QUICKBIRD, Remote 
Sensing, Urban Cartography, Urban Monitoring, 
Sharpness, Spatial 


1. ANTECEDENTES 


Ante la necesidad de contar con un sistema de 
información municipal que sustentara el cobro 
justo y oportuno del impuesto predial en la ciudad 
de Mexicali, surge el proyecto "Actualización del 
Registro Inmobiliario Municipal” en el marco del 


convenio de colaboración entre el XVI Ayunta- 
miento y la Universidad Autónoma de Baja Califor- 
nia (UABC). En este convenio, la UABC adquirió la 
responsabilidad de diseñar e implementar un sis- 
tema basado en el predio como unidad espacial de 
análisis básico. 

El predio es un limite espacial dinámico, es decir, 
en el tiempo, cambia de forma y dimensión, sea 
por fusión o subdivisión, o bien, cambia de conteni- 
do, sea por demolición, ampliación o nueva cons- 
trucción. La velocidad de estos cambios, en el 
caso de Mexicali, rebasaron la capacidad de la ofi- 
cina de Catastro de mantener actualizada la carto- 
grafía. 

El predio es también un límite subjetivo que 
demarca una propiedad, y que en muchos casos, 
como sucede en los asentamientos irregulares o 
suburbanos, aún cuando el predio existe física- 
mente, este no se encuentra documentado oficial 
mente, o existe discrepancia con lo registrado. 
Esto resulta un sobre o sub registro de propieda- 
des y por lo tanto en la asignación y cobro injusto 
del impuesto predial, 

Por lo anterior, en materia de cartografía se persi- 
guieron tres principales metas: la integración del 
total de predios urbanos de Mexicali (aproximada- 
mente 200,000) con precisión y rapidez: la identifi- 
cación de predios con incongruencia entre registro 
oficial documental registro espacial real, y final- 
mente, el diseño de un instrumento que permitiera a 
Catastro llevar a cabo un proceso de actualización 
permanente de la cartografía predial integrada por 
la Universidad. 


2. INTRODUCCIÓN 


En la elaboración y actualización de cartografía 
urbana, la imagen de satélite es un instrumento de 
gran utilidad, ello debido a las ventajas que ofrece 
su utilización con respecto a otras fuentes de 
información espacial. Algunas ventajas son las 
siguientes: la imagen satelital proporciona una 
visión completa y homogénea de la ciudad; al ser 
adquirida en secciones de tiempo relativamente 
cortos proporciona información relevante para el 
monitoreo urbano; el manejo informático de los 
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datos abre la posibilidad de realizar operaciones 
matemáticas y estadísticas complejas. 

Para obtener precisión y calidad en la elaboración 
de cartografía predial urbana es necesario utilizar 
como referencia una imagen con la mayor definición 
posible, es decir, una imagen donde los limites de 
los elementos urbanos son detectados con facili- 


Las técnicas de fusión, pretenden integrar gru- 
pos de datos con diferentes caracteristicas espa- 
ciales y espectrales para obtener productos de 
mejor calidad que los originales. (Vergara, 1997). 
De esta manera, al involucrar diversas fuentes de 
datos, es posible llevar a la información 
más allá de sus propias limitantes 
temporales, espaciales y espectrales. 

Existen diversas técnicas de fusión de datos: pero 
de manera podemos clasificarias de acuer- 
do a la complejidad del procedimiento que utiliza en 
básicas y complejas. Las técnicas básicas, se 
basan en procedimientos simples como las sustituti- 
vas (componente principal y canal de intensidad) y 
las aritméticas como la multiplicativa. Las técnicas 
complejas, se sustentan en combinaciones aritmét- 
cas como Broovey Transform, PBIM y SFIM (Liu, 
2000), o como las basadas en la teoría de Demp- 
ster-Shafer (Le Hégarat-Mascle et al 2003) y el 
modelo de Markov (Ersboll el al, 1998), algunas son 
onginadas a partir de múltiples descomposiciones y 
refleraciones como la Wavelet (Zhou et al, 1998), 
entre otras mas. 

De acuerdo a la revisión realizada por Zhang 
(2002) es común que la imagen obtenida por la apli- 
cación de las técnicas presente problemas de dis- 
torsiones del color y debilidad en los tonos, entre 
otros, asociados principalmente a la amplitud del 
rango espectral de la imagen pancromática, a la 
fusión de imágenes con resoluciones extremada- 
mente diferentes y de fechas distintas. Cabe pre- 
guntarse ¿la simple aplicación de las técnicas no 
garantiza la obtención de una imagen adecuada?, 
o bien ¿imágenes de un mismo sensor con proxi- 
midad de resolución y similar fecha de adquisición 
producen una imagen adecuada independiente- 
mente de la técnica aplicada? 

Con la intención de responder a los anteriores 
cuestionamientos y para seleccionar la técnica mas 
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conveniente para la elaboración de la cartografía 
predial de Mexicali, el objetivo del presente trabajo 
es aplicar y comparar tres técnicas de fusión de 
datos básicas (multiplicativa, sustitución del com- 
ponente principal, sustitución del canal de luminosi- 
dad) en una porción significativa de la imagen sate- 
iital de la ciudad de Mexicali. 

Para lograr lo anterior, se fusionan dos imágenes 
adquindas por un solo sensor, el QuickBird, en la 
misma fecha (Tabla 1) y de resolución próxima: a) 
imagen multespectral integrada por las bandas 
azul-verde-roja (BGR) resolución espacial de 240 
m, y b) imagen pancromática con resolución de 0.60 
m. Del proceso de fusión, se obtienen imágenes nue- 
vas, las cuales se examinan, y se caíífica la calidad 
visual que presentan, es dear, la clandad en el des- 
pliegue de los elementos urbanos en distintos usos 
de suelo, en diversos niveles de acercamiento. 

El software utilizado en el procesamiento de imá- 
genes es Idrisi32, desarrollado por el laboratorio 
de la Universidad de Clark, en Estados Unidos de 
Norteamérica. 


3. METODOLOGÍA 


Para cumplir con el objetivo del presente trabajo, 
se realizaron las siguientes tareas: el procesa- 
miento preliminar de las imágenes originales, la 
aplicación de las técnicas de fusión de datos y 
finalmente, la calificación y comparación de las 
imágenes resultantes, 


3.1. Procesamiento preliminar 

Inicialmente se realizó la homogeneización de 
las imágenes en cuanto a su referencia geográfi- 
ca, es decir, se registran en la proyección WGS84, 
zona 11 norte. Aunque la ciudad de Mexicali es 
topográficamente plana, se corrigió la imagen 
tomando como referencia las redes de puntos geo- 
désicos municipal y estatal existentes. 

De la cobertura total de la imagen, se extrajo el 
área de experimentación de aproximadamente 
1471 has, la cual se caracteriza por contener los 
usos de suelo de interés. Esta porción se encuentra 
delimitada por las coordenadas: 

Superior izquierda: 641051.40, 3616369.80; 
Inferior derecha: 644889 60, 3612537.60 

La Tabla 2 muestra los nombres asignados a los 
archivos. 


MN o 


Tabla 1. Características de las imágenes originales 
a PS A ATAR TS 
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12 ital Globe _| del 2003 

A Y — os 
450-900 Digital Globe_| del 2003 








3.2. Aplicación de las técnicas 

de fusión de datos 

Se obtuvieron imágenes secundarias como pro- 
ducto de la aplicación de las técnicas de fusión de 
datos que a continuación se describen, 


a. Técnica multiplicativa 

Consiste en la sobreposición multiplicativa de 
las imágenes originales (Esquema 1). Esto signifi- 
ca que cada número digital (ND) que compone la 
imagen pancromática es multiplicado por el 
correspondiente ND de cada banda de la imagen 
multiespectral, 


Esquema 1. Técnica multiplicativa 





De esta manera, se obtienen tres imágenes 
secundarias (PMS1, PMS2 y PMS3) cuyos valores 
espectrales se encuentran magnificados, de tal 
forma que el intervalo de valores originales de O a 
255 se amplificó de O a 65025. Las medias aritmé- 
ticas de los ND, se incrementaron como se mues- 
tra en la Tabla 3, 

Para regresar la imagen a valores aproximados 
a los originales, se lleva a cabo un ajuste de valo- 
res espectrales, que consiste en dividir los datos 
entre una constante, 


b. Técnica de sustitución del canal 

de luminosidad (intensidad) 

Esta técnica consiste en transformar las bandas 
originales del espacio de color BGR (azul-verde- 
rojo), al espacio de color constituido por los cana- 
les LHS (luminosidad-matiz-saturación). Poste- 
rormente se sustituye el canal de luminosidad (L) 
por la'imagen pancromática previamente ajusta- 
da; y se ejecuta el proceso inverso, es decir, se 
transforma el espacio LHS a BGR. (Vergara el al, 
1997). 
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Esquema 2. Técnica de sustitución 
del canal de luminosidad 


B (PHS1) 
/ 


Este proceso genera nuevas bandas (Esquema 2) 
cuyos valores resultan mayores a los onginales. Las 
medias de las imágenes secundarias son los 
siguientes B(PHS1)=132, G(PHS2)=132 y 
R(PHS)=136; por lo que es necesaño ajustar los 
valores, como en la técnica anterior, para aproximar- 
los a los originales indicados en la Tabla 3. 


G (M52) 


R (MS3) 





c. Técnica de sustitución 

del componente principal 
El análisis de componentes principales (ACP) es 
una técnica de transformación lineal. Dada una 
serie de bandas de imagen, el algoritmo ACP 
genera un nuevo grupo de imágenes, conocidos 
como componentes, no correlacionados entre sí, 
ordenados en términos de la cantidad de varianza 
que explican del conjunto de imágenes originales 
(Eastman, 1999:41), 

El primer componente contiene la mayor canti- 
dad de varianza de los valores de reflectancia cap- 
tados por el sensor, y los componentes posteriores 
se integran principalmente por ruido o por informa- 
ción no relevante, 

El algoritmo ACP, se aplica a la imagen multies- 
pectral para obtener tres componentes principa- 
les: C1, C2 y C3. La imagen pancromática se ajus- 
ta alos valores espectrales del primer componente 
(C1) y se ejecuta el proceso inverso de ACP para 
regresar los datos a su lectura onginal, sustituyen- 
do el primer componente por la imagen pancromá- 
tica (Chávez, 1991), 


Esquema 3. Técnica de sustitución 
del componente principal 


B (PC1) 
G (PC2) / k 
R (PC3) 


B (MS1) 
G (MS2) 





Una vez aplicado el algoritmo ACP. se obtuvo 
que el primer componente, como puede 
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Tabla 4. Porcentajes de varianza y valores de los componentes princi ae 
| COMPONENTE | C1_ | 







observarse en la Tabla 4, explica el 98% de la  quela Figura 1. Se observa una imagen mas o0scu- 

vanación de los datos. ra que la ongmal, y es posible percibir una mejor 
Los valores de las imágenes resultantes definición de los elementos urbanos 

B(PC1), G(PC2) y RIPC3), presentan promedios 

de 77, 138 y 122 respectivamente, por lo que es 

necesano ajustarios aproximadamente a los valo- 

res onginales especificados en la Tabla 3. 


3.3 Identificación de elementos urbanos 

Con las imágenes secundarias producto de apli- 
car las técnicas anteriores, se crearon composicio- 
nes en color verdadero. 

De las composiciones resultantes se tomaron 
secciones muestra de los usos de suelo de inte- 
rés (agricola, equipamiento y servicios, industrial, 
habitacional densidad baja y media), sobre tas 
cuales se realizaron desplegues a diferentes nive- 
les de acercamiento (zoom). En cada nivel se call- 
fica la definición de cada elemento urbano (viall- 
dades, manzanas, vegetación, predios y cons- 





trucciones) Figura 2. Composición técnica multiplicativa 
En la Figura 1 se muestra un acercamiento de la 
zona de uso habitacional densidad baja, en la La Figura 3 muestra la composición elaborada a 


composición en color verdadero elaborada a par- partir de las bandas producto de la técnica de sus- 
tir de las bandas originales. En este nivel de acer. — titución del canal de luminosidad. Es posible apre- 
camiento, los límites de los elementos son difu-  ciarleves deficiencias en la definición de la vegeta- 
$05 ción, sin embargo, la resolución espacial! se incre- 
mentó considerablemente con respecto a la técni- 
ca mulbipbcaliva 





' amposición con bandas originalos 





En la Figura 2 se muestra la composición elabo- 
rada a partir de las bandas producto de la técnica Figura 3, Composición técnica sustitución 
multiplicativa, con similar nivel de acercamiento del canal de luminosidad 
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En la Figura 4 se muestra la composición pro- 
ducto de la sustitución del componente principal 
por la imagen pancromática. Esta imagen presen- 
tó, al igual que la anterior, un incremento significa- 
tivo en la resolución espacial. Los elementos urba- 
nos se definen con claridad y no presenta proble- 
mas en las tonalidades de los objetos 





mM + 


Figura 4. Composición técnica sustitución 
del componente principal 


4, RESULTADOS 
»t 


La visualización de los elementos urbanos en los 
distintos usos de suelo, permitió valorar la calidad 
de la imagen resultante, Aún cuando las lécnicas 
tuveron como resultado un ennquecimento de las 
imágenes originales, este se manifestó en grados 
distintos 

Se observó una mejoría leve como resultado de 
la técnica multiplicativa, En cambio, las técnicas 
de sustitución (componente principal y luminosi- 
dad) aumentaron considerablemente la calidad 
visual de las imágenes originales, ello se percibe 
aún en niveles altos de acercamiento, donde se 
observan con claridad los limites de los elemen- 
tos. 

Entre las técnicas sustitutivas, la de huminosi- 
dad, presentó deficiencias en la definición de 
vegetación y vialidades no pavimentadas, esta 
imagen puede sin embargo, ser enriquecida con 
procesamientos posteriores para mejorar su des- 
pliegue. 

Con la técnica de sustitución del componente 
principal se obtuvo la mejor cahdad visual de la 
información. Las imágenes, producto de la aplica- 
ción de esta técnica, presentaron gran enriqueci- 
miento espacial de los datos, sin distorsionar los 
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valores, ni las tonalidades de los objetos, por ello, 
para su utilización en la elaboración de cartografía, 
aún cuando se genera a partir de una técnica bási- 
ca, no es necesario realizar posteriores procesa- 
mientos. 

La técnica de sustitución del componente princi» 
pal, se aplicó a la cobertura total de la imagen, 
(área urbana de Menucali) y a partir de ella, se ela- 
boró cartografía predial urbana en formato digital 
En la Figura 5 se observa el trazo de la cartografía 
de predios y construcciones sobre una porción de 
la imagen 





Figura 5. Cartografía predial urbana sobrepuesta 
a la imagen satolital 


5. CONCLUSIONES 


Aunque la fusión de datos permitió Hevar la 
información original mas allá de sus limitantes 
espaciales y espectrales, las técnicas aplicadas 
presentaron diferencias en la definición de ele- 
mentos urbanos, por lo que, algunas de ellas 
requieren de procesamientos previos o posterio- 
res. 

De esta forma, es posible comprobar que la 
homogeneidad del sensor y el manejo de imáge- 
nes de resoluciones y fechas próximas, en el pre- 
sente caso, no influyeron de manera signihicativa 
en evitar la distorsión en el color y la pérdida de miti- 
dez, en la imagen resultante. Por esta razón consi- 
deramos que el resultado obtenido se relaciona 
principalmente con la técnica aplicada. 

Las técnicas básicas son relativamente sencillas 
de aplicar, y aun cuando se trata de imágenes con 
rango espectral ampho, como es el caso de las imá- 
genes QuickBird, producen excelentes imágenes 
secundarias, e incluso en menor tempo que las 
técnicas complejas. 
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Para procesar las imágenes del proyecto de 
registro inmobiliario, se seleccionó la técnica de 
sustitución del componente principal, la cual per- 
mitió obtener la imagen de Mexicali con resolución 
de 0.60 m en color verdadero, como un excelente 
insumo en la elaboración de cartografía predial 
urbana y el monitoreo urbano. 

Finalmente, en la oficina de Catastro, se instauró 
un proceso de actualización permanente y norma- 
lizado de cartografía predial para alimentar al sis- 
tema, sustentado en la incorporación de imágenes 
satelitales fusionadas adquiridas en periodos 
específicos y en la verificación del crecimiento 
urbano como rutina de trabajo, para garantizar en 
el tiempo, un cobro más justo del impuesto predial. 


Catastral con 
zación de imagen Satelital en Menicali, B. C, ind 
XI Reunión Nacional SELPER México. 


* Zhang (2002) muestra la diferencia en el rango espectral 
diferentes sensores, 


ción de las imágenes se utilizó el total de puntos de ambas 
redes. considerando la coordenada del punto fis camente posi- 
cionado contra la coordenada del punto visualizado en la ima- 


gen. 
* En ldrisi32. a diferencia de otros 

indispensable 

elaborar 


programas (software), es 
realizar el ajuste de datos al rango 0 255 para 
en color. 
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Estimación de la diferencia entre datos de temperatura 
superficial del mar de imágenes AVHRR y mediciones in situ, 
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El objetivo fue examinar datos de temperatura 
superficial del mar en la Bahía de Campeche, oble- 
nidos de dos fuentes: i) datos medidos in situ, 
durante un crucero oceanográfico del 8 al 15 de 
febrero de 1996, con base en una red de 27 esta- 
ciones ubicadas mente, cubrnendo una 
área aproximada de 34,000 km2,; li) imágenes del 
Advance Very High Resolution Radiometer del 
satélite National Oceanic and Atmospheric Admi- 
nistration cercanas al perñodo del crucero. Las 
correlaciones entre las mediciones in situ y la infor- 
mación generada con las imágenes no fueron 
mayores que 0.75. Sin embargo, los patrones 
resultantes de distribución de temperatura superfi- 
cial del mar de las dos fuentes fueron similares cua- 
htativamente. Los cortes verticales refuerzan el 
comportamiento superficial de temperatura en 
ambos casos. Se detectan algunos errores imputa- 
bles al efecto atmosférico. Los resultados reafir- 
man que las imágenes son una manera práctica 
para estimar el campo de distribución de tempera- 
tura superficial del mar, semejante al que se obtie- 
ne con mediciones in situ. 


Palabras Clave: AVHRR, Bahía de Campeche, 
in situ, Isotermas 


ABSTRACT 


The objective was to examine sea surface lem- 
perature data in the Campeche Bay obtained from 
two sources: ¡) data measured ín situ during a 
research cruise carry out from February 8” to 15”, 
1996; the track was based on a grid of 27 stations 
located geographically and covering an approxi- 
mate area of 34, 000 km”; li) images of the Advance 
Very High Resolution Radiometer on National 
Oceanic and Almospheric Administration satellite 


near to the period of the research cruise. The 
correlations between the measurements in situ 
and the information ted with the images 
were not greater than 0.75. However, the patters of 
sea surface temperature distribution from the two 
sources were qualitatively similar. Vertical sections 
support the behavior of surface temperature in 
both cases. Some errors are detected attributable 
to the atmospheric effect. The results reaffirm that 
images are a practical way to estimate the field of 
distribution of sea surface temperatures similar to 
the estimate obtained by measurements in situ. 


Keywords: AVHRR, Bahía de Campeches, in 
situ, isotherms 


1. INTRODUCCIÓN 


Las | oceanográficas en México se 
incrementaron a partir de 1981 con la adquisición 
de los buques oceanográficos “Justo Sierra” y “El 
Puma” por parte de la Universidad Nacional Autó- 
noma de México (UNAM). La facilidad de contar 
con otra herramienta complementaria desde el 
aspecto de la oceanografía fisica, surgió con la ins- 
talación en 1994 en el Instituto de Geografía de la 
UNAM de una antena receptora de imágenes, 
Advance High Resolution  Radiometer 
(AVHRR) del satélite National Oceanic and 
Almosphenc Administration (NOAA), lo cual hizo 
posible comparar los patrones de distribución de 
temperatura superficial del mar (TSM), de una 
misma área geográfica derivados de dos fuentes. 

Comúnmente las mediciones de la TSM in situ se 
realizan con base en una red de estaciones, con una 
limitante, las mediciones no corresponden a un 
mismo punto en el tempo, entre mas grande es la 
red mayor es la diferencia en el tempo. Las imáge- 
nes no benen este mconveniente pero la exactitud 
del valor de TSM es menor, este problema 
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ha sido largamente discutido desde los años 1970 's 
(Li et al., 2001) No obstante las imágenes AVHRR 
son una herramienta de gran utilidad, por medio de 
ellas se logra observar la estructura de TSM con 
extenso alcance espacial, donde se pueden detectar 
fenómenos fisicos tales como plumas térmicas (Tang 


námica regional (Gibbs y Shaw, 2002), su restricción 
es que en ngor sólo son representativas de la tempe- 
ratura en una capa extremadamente delgada, sin 
embargo, permilen apreciar la distribución en la 
columna de agua de otros parámetros relacionados 
con la temperatura, siempre y cuando se verifique 
una buena correlación entre los valores de tempera- 
tura estimados a partir de los sensores remotos y los 
valores oblenidos in situ (Perez-Marrero et al., 1999). 
En este trabajo se aplicaron diferentes criterios 
para oblener compuestos de TSM a partir de un 
conjunto de imágenes, se calculó la correlación 
entre los datos de temperatura derivados de estos 
compuestos y los datos medidos in situ para esti- 
mar cual criterio proporciona mayor exactitud. Se 
hicieron cortes horizontales de TSM con los datos 
medidos in situ utilizando diferentes métodos de 
interpolación. Los compuestos y los cortes hor» 
zontales se compararon cualitativamente para 
determinar cuales coincidian más. Se cotejaron los 
valores de temperatura satelitales y los medidos in 
situ con mayor correspondencia en tempo, con la 
finalidad de analizar su exactitud. Se graficaron cor- 
tes verticales para observar su relación con los 
patrones de distribución de TSM. 
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2. MATERIALES Y MÉTODOS 


Mediciones in situ 

El área de estudio está situada al sur de la Bahía 
de Campeche, la cual es una de las regiones más 
preciada del mar patrimonial mexicano por el valor 
económico de sus pesquerías y la existencia 
de mayor importancia de depósitos de hidrocarbu- 
ros, recursos estratégicos en la economía nacio- 
nal y mundial. Por estas razones en la Bahía de 
Campeche se han desarrollado numerosos pro- 
yectos de investigación con carácter multidiscipli- 
nario, entre elos proyectos de oceanografía fisi- 
ca, en especial los relacionados con la hidrodiná- 
mica debido a sus implicaciones sobre las pes- 
querías, la distribución de desperdicios industria- 
les y la construcción de plataformas petroleras. 
Dentro de este marco se realizó el proyecto "Moni- 
toreo de prerreciutas de especies estuarino 
dependientes, comunidades bénticas y mareas 
internas del sur del Golfo de México” (MOPEED). 
En el ámbito ce este proyecto se efectuó el cruce- 
ro oceanográfico MOPEED 17, durante el cual se 
hicieron 27 estaciones oceanográficas distribui- 
das en 5 transectos perpendiculares a la costa (Fi- 
gura 1). Los transectos se numeran del uno al 
cinco (de oeste a este). El área de estudio se loca- 
liza entre 180 20" y 190 50' de latitud norte y 910 
30' y 940 30' de longitud oeste, con una superficie 
aproximada de 34,000 km2. 

La ubicación espacio temporal de las estaciones 
se realizó con el sistema de satélites Global Posi- 
tioning System (GPS), mediante el navegador de 
satélite TRANSITIGPS, con una precisión de 1 
segundo en el tiempo y 37 m en el espaco (Mag- 
navox, 1987). Se hicieron observaciones en 
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las 27 estaciones de la red con una sonda que mide 
conductividad, temperatura y presión (CTD) con 
exactitud de 0.001 mS/cm, O. 001 “€ y 0.05% FS res- 
pectivamente (www. generaloceanic.com, 2003), 
Las mediciones se iniciaron unos segundos antes 
de sumergir el CTD en el agua con el propósito de 
oblener el dato de TSM. El intervalo de muestreo 
fue de un décimo de metro en la columna de agua. 
Mediante el paquete computacional (paquete) del 
CTD (General Oceanics, 1990) se calculó salinidad 
y profundidad. 


Planos horizontales de temperatura 

de datos in situ 

Las isotermas de los cortes horizontales super 
ciales se trazaron mediante tres métodos de inter- 
polación. Á continuación se describen los proce- 
dimientos. 








ESTIMACIÓN DE LA DIFERENCIA ENTRE DATOS 


E 
ANI) FO 


13 


(1) Trazo manual combinado con los paqueles 
ARCIInfo (1992) y GRASS (1989). Se localizaron 
las estaciones sobre un plano de la zona, con los 
valores correspondientes de TSM para cada esta- 
ción se calcularon valores interpolados linealmen- 
te y se trazaron manualmente las isotermas, a par- 
tir de esta gráfica con los paquetes ARC/info y 
GRASS se obtuvo la disitnmbución horizontal de 
TSM que se muestra en la Figura 2(a), 

(11) Procesamiento utilizando los paquetes 
ARCAnfo y GRASS. Conlos datos de posición geo- 
gráfica, TSM y el paquete ARC/Info se graficaron 
las estaciones y se trazaron las isolineas, este 
paquete interpoló con el método inverso del cua- 
drado de la distancia. Finalmente con GRASS se 
elaboró la gráfica que se muestra en la Figura 2(b). 

(im) Procesaruento con el sistema SIGO (Padi- 
lhla-Pilotze, 2002). Con este sistema que 


23,1 


nn: ns mí sl AN 


(b) 


Figura 2. Distribución de TSM, datos del CTD, procesamiento (a) manual combinado con los paquetes 
ARC/nfo y GRASS y (b) con paquetes ARC/Info y GRASS. 
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utilizó el método de interpolación bilineal se obtuvo 
la distribución horizontal y vertical de temperatura. 
La Figura 3(a, b) corresponde a la distribución hor» 
zontal de TSM y temperatura a 10 m de profundi- 
dad. La Figura 4(a, b) presenta la distribución verti- 
cal de temperatura en los transectos | y Il. 


imágenes AVHRR 

EJ AVHRR está montado en los satélites de órbita 
polar NOAA y tiene una resolución de 1,1 km2 
sobre el nadir. Los datos de radiancia y brillantez de 
TSM derivan de los 5 canales con los que cuenta, 
dos en el visible (0.58 -0.68 m y0.72-1,10m) y 
tres en el infrarrojo (3.55 - 393 m, 10.3 -113m y 
115-125 m) Para los canales del infrarrojo el ruido 
aproximado es de 0.12*C (SeaSpace, 1995). 

Las imágenes AVHRR se capturaron por medio 
de la antena receptora instalada en el Instituto de 
Geografía de la UNAM. Éstas tienen una alta reso- 
lución radiométrica dentro de un intervalo de O 
a 1023(10 bits = 210). 

Los datos crudos que se oblienen de las imáge- 
nes se procesaron con las funciones indicadas en 
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el manual del sistema TeraScan (SeaSpace, 
1995), realizando el siguiente procedimiento: se 
extrajeron los datos AVHRR de los datos High 
Resolution Picture Transmission (HRTP), se cali- 
braron los datos y se convirtieron a temperatura 
radiométrica, en grados Celsius, se corrigieron los 
datos erróneos y se llenaron los huecos por inter- 
polación, se redujo el ruido y se calculó la TSM con 
el algoritmo MCSST (McClain et al., 1885). Poste- 
normente, con el sistema GRASS (1989), las imá- 
genes se corrigieron geométricamente y se 
remuestrearon con celdas de a 1 x 1 Km. 


Selección de | 

De las 12 imágenes AVHRR existentes cercanas 
en el tiempo al periodo durante el cual se realiza- 
ron las estaciones oceanográficas, se selecciona- 
ron ocho que tenian menor cobertura de nubes. 
La localización en el tempo tanto de las estacio- 
nes (lineas verticales) como de las imágenes 
AVHRR (círculos) se ve en la Figura 5. 

La Figura 6 muestra un contador de la frecuencia 
de datos (de O a 8) por pixel de las ocho imágenes. 





(b) 
Figura 3, Distribución horizontal de tamperaura, datos del CTD procesados con el SIGO. 
(a) superficial y (b] a 10 m de profundidad. 
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Figura 4. Distribución vertical de termperatura, datos del CTD procesados con al S1GO. 
fa) en el transecto | y (bi en el transecto ll, 
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Figura 5. imágenes y estaciones ubicadas en el tiempo. 
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Compuestos 

A partir del juego de las ocho imágenes con dife- 
rente cobertura de nubes, donde cada pixel puede 
tener de O a 8 datos, se decidió hacer tres com- 
puestos utilizando uno de los tres crtenos siguien- 
tes para cada uno de ellos. 

(1) Compuesto de TSM máximo (CM). En cada 

pixel se seleccionó el máximo valor de TSM, 
(1) Compuesto de TSM promedio (CP). Se calculó 
el promedio de los valores de TSM en cada píxel. 
(1) Compuesto de TSM tiempo (CT). Se tomó una 
imagen base por fecha y se fueron llenando los 
huecos correspondientes a la presencias de 
nubes con datos de las otras imágenes dando prio- 
idad a las imágenes más cercanas en tiempo a las 
fechas y horas en que se realizaron las estaciones 
oceanográficas. 

En la Figura 7(a, b, c) se presentan los compues- 
tos correspondientes. 


3. RESULTADOS Y DISCUSIÓN 


Si se compara los patrones de distribución de 
TSM de los tres compuestos y los tres productos 
derivados de las mediciones in situ observamos 
que: 

La gráfica procesada con los paquetes 
ARC/Info y GRASS (Figura 2(b)) muestra 
núcleos de agua fría y de agua caliente, que no se 
detectan en los otros patrones de distribución de 
TSM., 

El CP (Figura 7(b)) expone una mayor distribu- 
ción de agua de menor temperatura (a la izquierda 


Figura 6. Datos de TSM por pixel, 


de la laguna de Términos) hacia mar adentro, que 
no se observa en los otros patrones de distribución 
de TSM. 

Se puede notar que cualitativamente los patro- 
nes de distribución de TSM que tienen más seme- 
janza entre sí son: el que se obtuvo con el procesa- 
miento manual combinado con los paquetes 
GRASS y ARC/Info (Figura 2(a)), el elaborado con 
el SIGO (Figuras 3(3)), el CM (Figura 7(a)) y el CT 
(Figura 7(c)). 

Si se analiza los contes verticales de los transec- 
tos | y Il Figura 4(a, b) donde se localizan los gra- 
dientes más fuertes de temperatura, se observa 
que éstos también son visibles en los cortes hori- 
zontales. Al comparan la Figura Ya) y 4(a, b) se 
nota que la estructura de TSM se refleja hasta los 
5 m de profundidad, lo que está de acuerdo con los 
resultados presentados por Pérez el al. (1999) A 
los 10 m Figura 3 (b) la estructura se manifiesta 
diferente. 


Correlación 

Con el fin de calcular la correlación (Spiegel et 
al., 2001) entre los valores de TSM medidos in situ 
y los generados de las imágenes correspondientes 
a las estaciones oceanográficas, se hicieron cuatro 
señes de valores a partir de las mediciones con el 
CTD y los derivados de los compuestos máximo, 
promedio y tiempo. En los compuestos no se 
encontraron valores para las estaciones 21, 22, 23 
y 24 por lo que éstas no se consideraron en las 
señes (Tabla 1). Enla Figura 8 se exponen las grá- 
ficas de las seres. 
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Figura 7. Compuesto de TSM. (a) máximo, (b) promedio, (c) tiampo. 
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Tabla 1. Series de valores de TSM generadas con el CTD y los compuestos: 
máximo, promedio y tiempo, 
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Figura 8. Serios de valores de TSM procedentes del CTD y de los compuestos: 
máximo, promedio y tiempo. 


Los resultados de los cálculos de ta correlación 
entre la sere del CTD y cada una de las señes de los 
compuestos se muestran en la Tabla 2. 


Tabla 2. Correlación 









CTD — OM 0.75 
CTD - CP |0.70' 
CTD-CT [0.61 


La correlación más baja fue CTD CT, lo que 
sugiere que nuestro experimento se realizó ante 
un fenómeno fisico en el que el tiempo (en un lapso 
de días) no fue significativo. En algunos 


fenómenos fisicos el tiempo es importante debido 
a la dinámica de los procesos, sin er , si la 
estructura es estacionaria el tiempo no es 
relevante. Hay que resaltar que de los tres 
compuestos, cualitativamente el CP mostró ser el 
menos simil con respecto a los patrones de 
distribución de TSM derivados de las mediciones 
con el CTD, no obstante, su correlación fue mayor 
a la del CT. Para hacer el CP se calculó el 
promedio de los valores de TSM en cada píxel, a 
pesar de que en cada pixel los valores no 

al mismo punto en el tiempo ni el 
número de datos es el mismo. Ninguna de las tres 
correlaciones fue mayor que 0,75. Los resultados 
de las correlaciones en este trabajo se consideran 
aceptables. Algunos errores definitivamente se 
deben a efecto atmosférico (Ruttenbero, 1981), 
estos errores han sido ampliamente discutidos 


isa 


SELPER 


Kumar et al. (2003) señalan diferentes 
sobre este tema (Hagan 1989; 


investigaciones 
Wick et al. 1992; Walton et al. 1998). 


Pasos de satélite 

Se analizaron las coincidencias que pudiera 
haber entre las horas de los pasos de satélite de 
las imágenes y las horas en que fueron realizadas 
las estaciones oceanográficas, con el fin de com- 
parar los valores de TSM, Se encontró que la esta- 
ción 20 fue la más cercana al paso del satélite. Las 
estaciones 19, 25, 26 y 27 se aproximaron al paso 
del satélite (Tabla 3). De 27 valores de TSM corres- 
pondientes a las estaciones únicamente se pudie- 
ron comparar 5 contra los valores de las imágenes 
debido a que la sincronía entre la hora del paso de 
satélite y la hora en que se realizó cada estación 
oceanográfica es una condición deseable pero difi- 
cil de lograr. Por otro lado aunque se pudiera con- 
seguir esta sincronía, la cobertura de nubes en las 
imágenes es un fenómeno sobre el cual no se 
bene control alguno. La diferencia más grande 
entre los valores de TSM se encuentra en la mayor 
coincidencia en tiempo entre la hora que se realizó 
la estación 20 y el paso del satélite (Tabla 3). Se 
descarta que esta diferencia sea resultado de una 
medición errónea del CTD, considerando que los 
valores de TSM de sus estaciones vecinas no 
difieren de más de 0.3*C. Además los valores de 
temperatura por debajo de la superficie de la esta- 
ción 20 refuerzan el valor superficial (Tabla 4). 

Un análisis de la relación entre mediciones in 
situ de TSM y la derivada de imágenes AVHRR indi- 
can qué es posible obtener una exactitud de 0.5 
yo o mejor (Wellington et a., 2001; Wooster el al., 
2001). Lo que concuerda con las diferencias (Ta- 
bla 3) excepto en la estación 20 cuya diferencia 
excede 0.1*C al0.5*C señalado. 

Li et al. (2001) encontraron que la diferencia 
entre mediciones in situ de TSM y la derivada de 
imágenes AVHRR fue 0.4 *C durante el día y 0.2*C 
durante la noche. Lo que se ajusta con las diferen- 
cias (Tabla 3), si se considera que las estaciones 
25 y 26 fueron noctumas (3.35 hrs., 5.33 hrs.) y las 


Tabla 3. Diferencia de TSM. 


AAA a ale 

Estación | del satélite | Dehora | CTD |iImagen| TSM*C M * 

19 | 6627 | 68507 | 180 | 33 | 231 | 02 J44_1_23685 

aw [es2 | 60 | 01 |236]|230| 06 JH — 

25 | 81.47 | 8050 | 097 | 233 | 233 | 00 | 
zos | 230 | 229 


01 
O 11228 | 227 | 01 






83.47 | 8555 
Y 18531 1 8555 
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estaciones 19, 20 y 27 fueron diurnas (12.56 hrs., 
14.08 hrs. y 7.2 hrs.), las menores as se 
observan en las estaciones noctumas. La diferen- 
cia en la estación 20 excede con 0.2 al 0,4 *C indi- 
cado. 


4. CONCLUSIONES 


Si bien la correlación entre los valores de TSM 
generados a partir de las mediciones in situ y los den- 
vados de las imágenes no fue mayor que 0.75, los 
patrones de distribución de TSM emanados de las 
imágenes resultaron similares cualitativamente, a los 
que se obtuvieron de las mediciones in situ. La corre- 
lación del CP fue mayor a la del CT. La técnica utiliza- 
da en el CP, empleada en otros campos para proce- 
sar imágenes, es una buena opción para estimar la 
TSM, aunque hay que considerar que fue el com- 
puesto que presento agua de menor temperatura 
que no se observó en los otros patrones de distribu- 
ción de TSM. El criterio que se aplicó para elaborar el 
CM sería el más recomendado debido a que su 
correlación fue la más alta, a reserva de comprobar 
este resultado en futuros ensayos. 

Se identificó un fuerte gradiente de temperatura 
entre los transectos | y Il que se manifestó en todas 
las estructuras de TSM presentadas en este trabajo. 

El patrón de distribución de la TSM se refleja 
hasta los 5 m de profundidad. 

Los planos horizontales de TSM que tienen 
mayor similitud con los compuestos fueron los que 
se obtuvieron con el SIGO y con el proceso manual 
combinado con los paquetes ARC/Info y GRASS. 

De las cinco diferencia de la TSM cercanas al 
paso de satélite y las tomadas in situ cuatro estu- 
vieron dentro del rango 00C 0,2*C. 

Las imágenes no pueden sustituir las medicio- 
nes de TSM que se realizan en un crucero oceano» 
gráfico desde un punto de vista cuantitativo y están 
limitadas a proporcionar sólo valores de TSM pero 
tampoco el barco puede remplazar a las imágenes 
AVHRR debido a su alcance espacial para detec- 
tar fenómenos fisicos de gran escala. Las imáge- 
nes y los cruceros oceanográficos son 





VOL. 30 N" 1 
POMO 300 


20 


herramientas complementarias. 

Un factor aleatorio es la nubosidad que se con- 
vierte en una limitante para obtener información 
puntual. 
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RESUMEN 


La región de la laguna La Picasa, al sur de la pro- 
vincia Santa Fe, ha sufrido recurrentes inundacio- 
nes en los últimos años que provocaron un 
aumento espectacular de la superficie de este 
cuerpo de agua. Este fenómeno impactó negali- 
vamente sobre la población, el sistema productivo 
y las vías de comunicación. La alternativa de miti- 
gación propuesta es la construcción del denomi- 
nado "Canal Norte”, cuya finalidad es evacuar 
excedentes hidricos hacia el rio Paraná. Los obje- 
tivos de este trabajo fueron: cuantificar las varia- 
ciones temporales de las coberturas hídricas en la 
cuenca propia, e integrar la información analizada 
en un sistema de información geográfica a fin de 
aportar datos para la calibración del modelo hidro- 
lógico utilizado para estimar aportes al canal. Para 
el logro de estos objetivos se utilizaron imágenes 
satelitales Landsat 5TM, datos del SIG 250 (IGM) 
y curvas de nivel (Dirección de Geodesia Provin- 
cial). 


Palabras clave: mitigación de inundaciones, imá- 
genes satelitales, sistemas de información geo- 
gráfica. 


ABSTRACT 


In the last years, the region of La Picasa pond, in 
the south of Santa Fe Province, undergone recu- 
rrent floods that caused a very significant increase 
in the pond surface level, 

This produced a negative impact on the popula- 
tion, on the productive system and on both the rou- 
tes and ralroads; this impact has not still been 
totally overcome. 

The construction of the denominated North 
Channel was proposed as an alternative of miti- 
gation with the purpose of conveying water 
excesses from La Picasa towards the Parana 


River, thus depleting the pond. The objectives of 
this study were to assess the temporal variation 
of the water coverage of the diflerent water 
bodies in the pond catchment, to observe the 
runoff of precipitated water and to integrate the 
analyzed information in a geographical informa- 
tion system in order to supply data for the calibra- 
tion of a hydrological model used to estimate late- 
ral inflows to the channel, 

For the study, ft was used Landsat 5TM Satellite 
Images, data of GIS 250 and maps of the area with 
elevation contours (source: Dirección de Geodesia 
Provincial). 


Keywords: flood mitigation, geographical infor- 
mation systems. 


1. INTRODUCCIÓN 


La laguna La Picasa está ubicada al sur de la pro- 
vincia de Santa Fe, en el límite con la provincia de 
Buenos Aires, abarcando parcialmente los distri- 
tos de Aarón Castellanos y Diego de Alvear (Bac- 
caglio et al 2003). La laguna forma parte de una 
extensa cuenca que lleva su nombre, que involu- 
cra las provincias de Santa Fe, Córdoba y Buenos 
Aires con una superficie total de 550.000 ha (Mi- 
nisterio de Obras, Servicios Públicos y Vivienda 
del Gobierno de Santa Fe, 2001). 

A partir de septiembre de 1997 se registraron 
en la zona lluvias extraordinarias que práctica- 
mente duplicaron los valores medios anuales de 
precipitación, pasando de 900 mm a 1700 mm 
aproximadamente. Como consecuencia de esto 
la laguna sufrió un incremento en su superficie de 
7.841 ha (año 1994) a 31.500 ha (en 2003) (Bac- 
caglio et al,, 2003), mientras que los correspon- 
dientes niveles que alcanzó la laguna aumenta- 
ron de 98.7 m hasta superar los 105 m en la 
actualidad (Ministerio de Obras, Servicios Públi- 
cos y Vivienda del Gobierno de Santa Fe, 
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2001, Baccaglio et al., 2003). 

El fenómeno de la expansión de la laguna La 
Picasa se debe a varios factores. El incremento en 
las lluvias es el principal de ellos pero también se 
reconocen otros factores importantes. Uno de 
ellos es el cambio histórico en el clima de la zona 
que pasó de ser templado subhúmedo hasta 
semiárido a marcadamente templado húmedo a 
partir de la década del "70. Esto produjo que en las 
últimas décadas los excesos de agua sean supe- 
nores a las salidas por evapotranspiración y dre- 
naje superficial, y que la capa freática ascienda y 
añiore en las zonas más deprimidas del paisa- 
je.(Tchilinguirian et al., 2003). 

Otro factor importante es que las condiciones 
naturales del sistema fisico de la cuenca influyen 
en la persistencia del anegamiento, ya que es una 
cuenca endorreica sin red de drenaje jerarquiza- 
da, en la cual la salida del agua es lenta porque se 
produce predominantemente por evapotranspira- 
ción o infiltración desde las zonas deprimidas. Por 
último deben mencionarse los factores antrópicos 
que introducen modificaciones al medio, La cuen- 
ca de La Picasa sufrió degradación en sus suelos y 
disminución de la capacidad de almacenamiento 
por el crecimiento de la superficie agrícola en detri- 
mento de la ganadera y por una sobreexplotación 
del suelo con prácticas intensivas de dobles culti- 
vos anuales. Á esto se agrega el incremento en la 
infraestructura de caminos, vías férreas, rutas y 
canales que alteran el normal escurrimiento de las 
aguas (Ministerio de Obras, Servicios Públicos y 
Vivenda del Gobierno de Santa Fe, 2001). 

Por último existen en el área canalizaciones o 
cuneteos, generalmente precarios y sin proyec- 
to, realizados para mitigar en forma rápida una 
situación de inundación excepcional, adquirien- 
do el carácter de trasvasamiento semi- 
permanente. Estos son de escasa eficiencia 
cuando se han construido por los motivos 
expuestos anteriormente pero, se han observado 
en el campo obras hidráulicas que presentan un 
perfil más permanente con obras, de arte de cons- 
trucción adecuada, que merecen “a posteriori” un 
análisis particular porque probablemente modifi- 
quen en forma definitiva la conformación del sis- 
tema hídrico actual (Instituto Nacional del Agua, 
2004), 

Entre las consecuencias más importantes de 
esta situación de excesos hídricos se mencionan 
las pérdidas en el sistema productivo y su impe- 
nosa necesidad de cambio o cambios en el siste- 
ma productivo, pérdida de superficie de pasto- 
reo, salinización de suelos, pérdida de caminos 
rurales, salinización de pozos de agua rurales y 
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el anegamiento de pozos sépticos por ascenso 
E de la capa freática (Tchilinguirian et al., 
3). 

Como alternativa de mitigación la Subsecretaría 
de Recursos Hidricos de la Nación (SSRH) propu- 
so la construcción de un canal de descarga deno- 
minado Canal Norte desde la laguna hasta el río 
Paraná, que permita evacuar los excedentes de 
agua. El proyecto de diseño de este canal fue enco- 
mendado al Instituto Nacional del Agua (INA). De 
toda el área que atraviesa el canal se prestará 
especial importancia al conjunto de lagunas inter- 
conectadas entre sí llamadas La Barrancosa, La 
Larga, Tuerto Venado, La Infinita, Martín Garcia y 
Los Patos, ya que le otorgan al sistema hídrico una 
gran capacidad de almacenamiento y amortigua- 
ción de los excedentes hídricos. (Subsecretaria de 
Recursos Hidricos, INA, 2004). 

Este trabajo forma parte de las diversas tareas 
encaradas por el INA en la etapa de prediseño del 
proyecto de construcción del citado canal de des- 
carga y para la evaluación de su impacto económi- 
co enla región. 


2. OBJETIVOS 


1-Cuantificar y evaluar la variación de la cobertu- 
ra hídrica de los siguientes cuerpos de agua entre 
las fechas febrero, abril, junio, septiembre y 
noviembre del año 2001: 

1,1-Cuerpos de agua de la totalidad de la cuenca 
de la laguna La Picasa. 

1.2-Cuerpos de agua de las subcuencas La Pica- 
sa-Los Patos, Los Patos-Quimo, Quirno- Juncal y 
Las Encadenadas. 

1.3-Cuerpos de agua de las micro cuencas defi- 
nidas dentro de la subcuenca Los Patos-Quirno. 

1.4-Cuerpos de agua de las lagunas La Barran- 
cosa, La Larga, Tuerto Venado, La Infinita, Martin 
García y Los Patos. 

2-Integrar la información generada en un Siste- 
ma de Información Geográfica para aportar 
datos que permitan contrastar y ajustar los resul- 
tados obtenidos con el Modelo de Simulación 
Es Continuo (HEC-HMS) utilizado por 
el CRL. 


3. MATERIALES Y MÉTODOS 


3.1 Área de estudio 

El área de estudio fue definida por el Centro 
Regional Litoral (CRL) del INA en base al análisis 
de curvas de nivel cedidas por la Dirección de 
Geodesia de la provincia de Santa Fe, encon- 
trándose al este de la laguna La Picasa, al 
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norte del límite con la provincia de Buenos Áires, al 
sur de la localidad de Venado Tuerto y al veste de 
las nacientes del arroyo del Medio. Para su análi- 
sis la cuenca fue subdividida en cuatro subcuen- 
cas: La Picasa-Los Patos, Los Patos-Quirmo, Quir- 
no- Juncal y Las Encadenadas (Figura 1). 

La subcuenca La Picasa-Los Patos fue a su 
O y dos micro cuencas (Figu- 
ra2) 


3.2 Información utilizada 

El estudio se llevó a cabo utilizando imágenes 
salelitales Landsat TM (cedidas gratuitamente, 
para su uso en emergencias, por la Comisión 
Nacional de Actividades Espaciales, CONAE), 
path 227, row 084 de las fechas 12 de febrero, 1* 
de abril, 20 de junio, 8 de septiembre y 11 de 
noviembre de 2001, El año 2001 fue elegido por 
ser un año en el cual se registraron 1412 mm de llu- 
via sobre un estado de alta saturación anteceden- 
te. Las fechas de las imágenes fueron elegidas en 
función de la ausencia de nubes que permiten su 
correcta visualización y procesamiento. 


3.3 Metodología de 
El procesamiento de las imágenes se inició- 
transformando los contajes digitales (DN) de cada 








iaa 
Figura 2: Detalle de las microcuencas de la 
subcuenca La Picasa-Los Patos. 


banda de las imágenes TM a radiancia a panir de: 
los coeficientes de calibrado correspondientes, 
suministrados en cada imagen. La ecuación lineal 
aplicada para obtener la radiancia es: 
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Le [móor sr.cm ES —— a 1 


donde: 
k: número de banda (k: 1, 2,3.465) 
Liz :walor digital máximo (254) 
-  tadiancia correspondiente a DN" = DN... 
a 1.5 
L”,., radiancia correspondiente a DN'= 0, para k: 
1.5 


Posteriormente se transformaron los valores de 
radiancia (L.(5.)) en reflectividad (R, (?).) a nivel del 
sensor (reflectividad a tope de atmósfera). Para la 
transformación se aplicó la siguiente expresión: 

: | 


nd.L o, 
Ev y. COSO 


' 


| 

donde: 

L: reflectividad en la longitud de onda considera- 
da (adimensional, entre O y 1), 

d: distancia de la Tierra al Sol, en unidades astro- 
nómicas, 

E,: irradiancia exoatmosférica, 

6 : ángulo cenital 


Rot) > 


En segundo lugar se realizó la corrección atmos- 
a a EE. 
grama por valores mínimos (Chávez, 1988). Este 
A 

la atmósfera a través de un cálculo que se hace 
directamente sobre la imagen para la determina- 
ción de la radiancia medida por el sensor sobre 
áreas oscuras. Esas áreas oscuras en teoría debe- 
rían ser negras (=0 por ciento de reflectancia) pero 
a causa de la dispersión atmosférica los corres- 
pondientes pixeles presentan un nivel digital distin- 
to de cero. Bajo estas hipótesis sólo una dispersión 
simple puede ser removida. 

Operacionalmente, la selección del nivel digital 
minimo apropiado para la implementación de la 
corrección puede ser oblenido por el histograma 
de frecuencias de la imagen digital (Chávez, 1996). 

Este algoritmo apunta principalmente al análisis 
de las bandas del espectro visible donde usualmen- 
te hay un marcado aumento en el número de pixe- 
les con valor distimo de cero, es decir, en escala de 
grises. La determinación del valor digital minimo 
para una banda en particular se basa en la hipóte- 
sis que asume el método en el hecho de que existe 
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una alta probabilidad de que por lo menos algunos 
pixeles dentro de una imagen sean negros. Esas 
superficies oscuras pueden corresponder a (1) 
áreas sombreadas por nubes o efectos topográfi- 
cos (baja reflectancia en todas las bandas), (2) cuer- 
pos húmedos (baja reflectancia en las longitudes 
de onda roja e infrarrojas), (3) áreas con densa 
vegetación 


tro electromagnético) y (4) superficies con combi- 
naciones de algunos de los factores anteriores 
(Kaufman y Sendra, 1988). 

La siguiente ecuación define el nivel de salida en 
la posición ¡,j para la banda k, DN ,, como el nivel 
digital original DN .,, menos nivel digital mínimo 
para la misma banda k, DN ,.,.... 





Este método fue utilizado por la ventaja de 
requerir sólo información contenida en la propia 
imagen. Se ha optado por este método ya que no 
se dispone de perfiles atmosféricos para la zona. 

Finalmente se corrigió geométricamente la 
escena, llevando su eferencia al sistema 
Gauss-Kruger, faja 5 (datum WGS84) mediante 
el método del "Vecino más Próximo”, ya que es el 
que menos modificaciones ocasiona al valor de 
reflectancia del pixel original. Se diseñó un mode- 
lo de procesamiento de imágenes satelitales (mo- 
deler) conteniendo dos indices calculados a par- 
tir de la reflectividad en diferentes bandas: NDVI 
y un índice de colorísaturación, más la banda 
azul, arrojando como resultado una imagen de 
tres bandas para su respectiva clasificación. 
Para la selección de pixeles se utilizó un método 
de clasificación híbrido o mixto (no supervisado / 
supervisado). Tal como se indicó, este método no 
supervisado no implica conocimientos previos 
del área de estudio, de modo que la intervención 
del investigador se orienta hacia la interpretación 
de los resultados; por tanto, en esta estrategia se 
asume que los DN de la imagen forman una serie 
de grupos o conglomerados (clusters) de relativa 
nitidez según sea el caso. Se seleccionaron entre 
70 y 80 números de clusters para una adecuada 
representación de la información. Estos clusters 
oestos grupos de pixeles son equivalentes en tér- 
minos de su comportamiento espectral más bien 
homogéneo y, por tanto, definen clases temát- 
cas de interés homogéneos. Estas forman parte 
de las muestras de entrenamiento del método 
supervisado. Se propone el método mixto Mibr- 
do, ya que el método supervisado puede 
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ornato afimrreneo subjetivo A 

puede carecer de significado fisico real 
naar 1997). Se efectuó un mapa temáti- 
co, con dos clases temáticas, agua superficial y 
encharcamiento, sobre el cual se realizó un proce- 
so de enmascarar (recode) para seleccionar sola- 
mente la clase de agua superficial. 

Se generaron cinco mapas multitemporales 
fusionando los mapas temáticos de clase agua de 
dos y tres fechas: febrero-abril, junio-septiembre, 

septiembre-noviembre, febrero-abril-junio, junio- 

septliembre-noviembre para la cuenca total. Sobre 
la imagen final de cada fecha se efectuaron 37 cor- 
tes correspondientes a las subcuencas y micro 
cuencas y se cuantificaron todos los cuerpos de 
agua presentes en hectáreas . 

Se digitalizaron y cuantificaron las coberturas 
hídricas en hectáreas de las lagunas La Barranco- 
sa, La Larga, Tuerto Venado, La Infinita, Martín 
García y Los Patos, ubicadas en la subcuenca La 
Picasa-Los Patos. Se elaboró un mapa para cada 
fecha comparando la superficie de todas las lagu- 
nas y un mapa por cada laguna comparada a tra- 
vés de todas las fechas. 

Los datos obtenidos fueron integrados en un sis- 
tema de información geográfica para ser usados 
como información para calibrar el Modelo de Simu- 
lación Hidrológica Continuo (HEC-HMS). 

Para el procesamiento de las imágenes se utili- 
zó el programa ERDAS (ERDAS) y para la digitali- 
zación de las cuencas y cuerpos de agua el pro- 
grama ARCVIEW (ESRI), ambos en un ambiente 
Windows PC. 
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4. RESULTADOS OBTENIDOS 


Se elaboraron tres Tablas con los datos de super- 
ficies obtenidos y un total de 42 mapas para facilitar 
la visualización e imerpretación de los datos. 

Se obtuvieron cinco mapas multitemporales de 
la cuenca total de la clase temática agua de dos y 
tres fechas: febrero-abril, junio-septiembre, sep- 
tiembre-noviembre, febrero-abril-junio y  junio- 
septiembre-noviembre de 2001 (ver ejemplo en 
Figura 3). 

Se obtuvieron 37 mapas multitemporales de 
las mismas dos y tres fechas y de la clase temáti- 
ca agua, de las 4 subcuencas y de las 33 micro- 
cuencas (ver ejemplo en Figuras 4 y 5). 

Se obtuvieron cinco mapas de la clase temática 
agua de las siguientes fechas: febrero, abril, junio, 
septiembre y noviembre de 2001 con el detalle de 
las lagunas La Barrancosa, La Larga, Tuerto Vena- 
do, La Infinita, Martín García y Los Patos, ubicadas 
en la subcuenca La Picasa-Los Patos, asi como 1 
mapa multitemporal de las 5 fechas juntas para 
cada laguna, siendo un total de seis (ver ejemplo 
en Figura 6). 

Se presentan las Tablas con los resultados de los 
cálculos de superficies de los cuerpos de agua de la 
cuenca total y de las cuatro subcuencas (Tabla 1), 
de las treinta y dos micro cuencas (Tabla 2) y para 
las lagunas La Barrancosa, La Larga, Tuerto Vena- 
do, La Infinta, Martín García y Los Patos (Tabla 3) 
para las fechas analizadas (febrero, abril, junio, sep- 
bembre y noviembre del 2001). 





Figura 3: Mapa multitemporal de la cuenca total para febrero y PP de 2001. 
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Figura 5: Mapa multitemporal de la microcuenca 18 para febrero, jurio y noviembre de 2001 
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Figura 6: Mapa multitemporal de las lagunas seleccionadas para abril y noviembre de 2001. 


Tabla 1: Suparhficia de los cuarpos de agua de la cuenca total y las subeuancas medidas an hectáreas. 
"SD"; sin dato. 


Superficie de los cuerpos de agua de la cuenca total y las subcuencas 





(Valores en hectáreas) 


CUENCA LUJTAL, 


LA PICASA-LOS P41O05 


AIN PATOS AMARNO 
OLIKNO JUNCAL 
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FERRERO 


5.048 


303 


4.147 


1H 1! JUXIO SEPTIEMBRE  —XOVIEMBRE 















8.069 10,976 
7674 3 252 7.974 9.193 
2.308 2.108 2,090 yO 
14,282 13,644 11,451 14,859 
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Tabla 2: Sup. de los cuerpos de agua de las micro cuencas medidas en ha. S/D: sin dato. 


Superficie de los cuerpos de agua de las micro cuencas 


Valores mm 
FEHHEMIO ABHI1 Mn 5SEPFTILVMHKI MIVIFAIHIKI 
hectares. | 

















103 
41 
1049 
179 





Tabla 3: Superficie de los cuerpos de agua de las lagunas La Barrancosa, La Larga. 
Tuerto Venado, La Infinita, Martin Garcia y Los Patos medidas en hectáreas. 
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Aplicación del modelo 
hidrodinámico SIS Flow 


A modo de reproducir el comportamiento hidro- 
lógico - hidráulico del tramo laguna La Picasa - 
laguna Los Patos en el año hidrológico extremo 
2001, se implemento el modelo hidrodinámico 
ISIS Flow al sistema Laguna La Picasa Laguna 
Los Patos Laguna El Chañar, utilizando la informa- 
ción topográfica de perfiles transversales confor- 
me a obra del canal construido y del camino veci- 
nal a la salida de la laguna Los Patos, relevados en 
el año 1999/2000 por la Dirección Provincial de 
Obras Hidráulicas 

Se simuló el comportamiento hidrológico del 
canal en el año 2001 y se compararon los resulta- 
dos obtenidos con aforos realizados por la Subse- 
cretaria de Recursos Hidricos de la Nación, aguas 
abajo de la sección de la Laguna de Los Patos. Los 
resultados estimados de la modelación fueron 


menores alos observados. Debe agregarse al cau- 


dal calculado, el caudal ingresado en el tramo Los 
Patos Sección de Aforo. 

También se utilizó para contrastar los resultados 
obtenidos del sistema modelado, las áreas cubier- 
tas por inundación en cada laguna a partir de la 
información procesada de imágenes satelitales 
tomadas los días 12/02, 01/04, 20/06, 08/09 y 
11/11 del año 2001 (Tabla 4). 

Se aplicó el modelo hidrodinámico a dicho año 
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hidrológico extremo y se determinaron los niveles 
alcanzados en las fechas indicadas en las imáge- 
nes. Á partir de dichos niveles se calcularon las 
áreas de cada laguna, asociadas a los niveles y 
se compararon con los resultados obtenidos de 
las imágenes. Las áreas asociadas a los niveles 
de las diferentes lagunas fueron calculadas a par- 
tir de las curvas de nivel de las planchetas ¡GM 
que tienen una equidistancia de 2.50 m y de los 
perfiles medidos de los relevamientos topográfi- 
cos mencionados anteriormente. Estos perfiles 
corresponden aproximadamente a tres por lagu- 
na y su ubicación planimétrica fue ajustada en 
forma aproximada, dado que sus coordenadas 
planimétricas planas se encuentran en Campo 
Iinchauspe, distinta a la faja a la usada por el pro- 
yecto. 

A partir de esta calibración aproximada y compa- 
rando los resultados obtenidos de las imágenes 
satelitales, puede observarse en la Tabla siguiente 
que se obtiene una aproximación razonable para 
los diferentes niveles de las lagunas y en los perio- 
dos de estiaje. 


A continuación se muestra en la Figura 7 la varia- 
ción temporal da las áreas inundadas de cada lagu- 
na (Tuerto Venado, La infinita, Martín García y 
Los Patos) de acuerdo a los resultados observa- 
dos y calculados por el modelo, con datos del año 
2001. 


Tabla 4: Comparación entre datos observados y calculados en el modelo hidrodinámico ISIS Flow. 
* D/O: Dato observado en imagen satelital, DIC: dato calculado en el modelo. 
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5. CONCLUSIONES 


Los valores de superficie hídrica para la cuenca 
total aumentaron en el considerado desde 
22.403 ha hasta 37.608, (67,87%), registrándose 
un cambio en la tendencia en el mes de septiem- 
bre . 

En las cuatro subcuencas analizadas también 
se registraron aumentos finales, pero las tenden- 
cias del cambio y los valores porcentuales variaron 
entre ellas. Para la Picasa-Los Patos el cambio fue 
de un 83.76% y la tendencia fue siempre en 
aumento si bien no se contó con información para 
el mes de junio. La subcuenca Los Patos-Quirno 
presentó aumentos y disminuciones alternativas 
entre los meses analizados y el porcentaje de 
aumento final fue del orden del 54.64%. Quirno 
Juncal y Las Encadenadas mostraron una tenden- 
cia similar que registró disminuciones en los 
meses de junio y sept y valores porcentua- 
les de aumento de 87.18% y 63.16% respectiva- 
mente. 

De las treinta y dos micro cuencas analizadas 
solo dos (6 y 7) registraron valores menores a los 
iniciales, 128 ha a 32 ha y 135 ha a 100 ha, respec- 
tivamente. 

De las treinta y un restantes, diecinueve de ellas 
presentaron una disminución de su valor de super- 
ficie en el mes de septiembre. Esto podría deberse 
posiblemente a como se distribuyeron las lluvias 
precipitadas en el área. 

Para las seis lagunas analizadas en particular 
también se registraron aumentos en sus valores 
finales excepto para Tuerto Venado que disminu- 
yó sensiblemente su superficie, pasando a ocupar 
de 489,33 ha a 486,99 ha. La Barrancosa aumen- 
tó su superficie un 23.33%, La Larga un 75.97%, 
La Infinita un 99.56% y Los Patos un 18.56%, 
registrando esta última una leve disminución de 
superficie para el mes de junio. 

La utilización de estos datos para el ajuste del 
modelo song utilizado resultó ampliamente 


El uso de información satelital permitió recupe- 
rar información que de otro modo no se hubiera 
obtenido ya que no se contaba con registros de 
campo de las coberturas hídricas para las fechas 
de interés. La elaboración de mapas multitempo- 
rales permite efectuar comparaciones rápidas y 
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pa y hacer un análisis de los cambios regis- 


La q A utilizada en el procesamiento 
de la información es de uso generalizado, no obs- 
tante, su integración y aplicación a un proyecto 
de diseño de una medida estructural, como lo es 
la obra de canalización de la descarga de la lagu- 
na La Picasa hacia el río Paraná es de gran 
actualidad y reviste pocos antecedentes en el 
país. 
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RESUMEN 


El presente trabajo plantea una metodología 
para evaluar la precisión de la ortorectificación de 
imágenes IKONOS (GEO, pan sharpened) de una 
zona urbana de Caracas. Esta metodología con- 
sistirá en ortocorregir la imagen original al Datum 
local de la Canoa (PSAD56), usando 50 puntos de 
control obtenidos de mediciones Diferenciales 
GPS (DGPS) y un DEM obtenido a partir de las cur- 
vas de nivel de las cartas topográficas a escala 
1:100.000 del Instituto Geográfico de Venezuela 
Simón Bolivar (IGVSB). 


SUMMARY 


The present paper presents a methodology to 
assess the ortho-rectification accuracy performed 
on IKONOS 1 meter pan-sharpened data collected 
over the north-western sector of Caracas city. This 
methodology consisted of re-projecting the original 
image to the local datum of La Canoa (PSAD 56), 
then, to apply a set of 50 ground control points 
(GCP"s) obtained by means of GPS differential 
observations (DGPS), as well as a 20 m DEM 
obtained from contour lines from existing topo- 
graphic chants at the scale of 1:100.000, produced 
by the Instituto Geográfico de Venezuela Simón 
Bolivar (IGVSB). 


1. INTRODUCCIÓN 


El objetivo de este trabajo consistió en plantear y 
evaluar la precisión de una metodología para la 
ortocorrección de imágenes del satélite IKONOS 
en una zona Urbana de Caracas. 

La metodología de ortorectificación de imágenes 
IKONOS que se ha utilizado en el Centro de Proce- 
samiento Digital (CPDI) del Instituto de Ingeniería 
(FI1) se basa en la elaboración de Modelos Digitales 
de Elevaciones (DEM) oblenidos de mapas topo- 
gráficos a escala 1:1.000 y 1:5.000, complementa» 
dos con puntos de control. Esto plantea serias res- 
tricciones para la ortorectificación de imágenes 
IKONOS en aquellas zonas donde no se disponga 


de cartografía a estas escalas. 

Por esto se consideró el estudio de otra metodo- 

| ao para ortocorregir imágenes de satélites de 
alta resolución espacial (1 m y 4 m). 

La metodología usada en el presente trabajo eva- 
Iwó ta precisión de la ortorectificación de imágenes 
IKONOS de im (GEO, pan sharpened) de una 
zona urbana de Caracas. Esta consistió en repro- 
yectar la imagen original al datum local de la 
Canoa (PSAD56), usando 50 puntos de control 
obtenidos de mediciones Diferenciales GPS 
(DGPS) y un DEM obtenido a partir de las curvas 
de nivel de las cartas topográficas a escala 
1:100.000 del IGVSB. 

Se procesó la imagen IKONOS usando el software 
OrthoEngine versión 7.0, este utiliza el modelo desa- 
rrollado por el Dr. Thierry Toutin del Canadian Centre 
lor Remote Sensing (CCRS), Natural Resources 
Canada (Toutin and Cheng, 2000). 

Para evaluar las precisiones de la técnica utili- 
zada, se desarrolló el estudio en una subescena 
de la imagen que comprende la zona Norte-QOeste 
y Norte-Centro de Caracas con un área de 1.170 
Hectáreas y con alturas entre 890 y 1.200 msnmm 
(Figura 5). También se evalúo la influencia de la 
cercanía de áreas montañosas, cuando se usa el 
modelo de ortocorrección. 

Los resultados se evaluaron correlacionándolos 

con los puntos de chequeo GPS, con los ortofoma- 
pas del Proyecto Cartocentro del IGVSB, así como 
con el producto ortocorregido con los mapas 
1: 1,000 y 1:5.000 en el CPDI (Ruiz, 2.000). 


2. ANTECEDENTES 


Atkinson y Smith (2.001) demostraron que los 
puntos de control terrestres (SCPs) originados de 
la campaña DGPS tienden a una mejor precisión 
en la rectificación de imágenes SPOT y TM que 
aquellos derivados de mapas topográficos (medi- 
dos por el RMS), concluyendo que la precisión de 
la rectificación depende de la precisión de los pun- 
tos de control GCPs. 

Gerlach (2.001) plantea que la elevación del sol 
y la inclinación del satélite sobre las diferentes 
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zonas de la escena tienden a crear variaciones en 
el contraste de la imagen cuando se trata de esce- 
nas pancromáticas sobre ciudades, porque se real- 
zan las sombras producidas por el sol y las zonas 
no vistas por el sensor, las cuales tienden a hacer 
que se pierdan detalles sobre el terreno, 

Main (2000) recomienda que la escogencia de 
los GCPs con alta calidad de reconocimiento en la 
imagen es de vital importancia para ubicar estos 
con alta precisión. Se sugiere la toma de puntos 
por ejemplo sobre la superficie asfáltica pintada, 
intersecciones de caminerias o en la “T” de una 
cancha de tenis. También se recomienda para 
obtener productos a escala 1:4.800 un RMS de 1,9 
pixeles. Además define un alto nivel de medicio- 
nes para la campaña GPS, en la cual, se deben 
obtener precisiones de 0,7 m en planimetria y 1,5 
m en altimetría. También recomienda que por cada 
10 Km de área se deben usar como mínimo 10 pun- 
tos de control GPS, 

Ganas, et al (2.001) realizó un trabajo en los sub- 
urbios de Atenas-Grecia que consistió en ortoco- 
rregir imágenes IKONOS GEO pancromáticas en 
un área plana, con variaciones entre los 70 y 250 
m, donde se usaron 17 puntos de control obteni- 
dos de una campaña GPS diferencial y un DEM a 
20 m obtenidos de los mapas topográficos de Gre- 
cia a escala 1:50.000, un error cuadrá- 
tico medio (RMS) de 0,6 pixeles y el mismo error 
para los puntos de verificación desde 15,97 pixe- 
lesa0, 9 pixeles. 

Dial (2.000) concluye con la experiencia de los 
mapas hechos desde fotografías aéreas que exis- 
te una constante en la resolución de la imagen, y la 
escala y precisión de los mapas. Se sugiere una 
precisión horizontal para ciertas escalas: Para 
mapas a escala 1:2.400 el radio de error en el 90% 
de la escena no debe ser mayor a 2,0 m (CE9O de 
2.0 m) y a escala 1:4.800 el CE9O debería ser de 
4,1 m. Para la primera, Corbley (2.000) sugiere 
una precisión de 3 m en altimetria para los puntos 
de control terrestre, 


IKONOS 

El satéite IKONOS puesto en órbita el 24 de sep- 
tiembre de 1.999 por la empresa Space Imaging Inc. 
Es el primer satélite comercial que ofrece imágenes 
de alta resolución espacial de 1 m en modo pancro- 
málico y 4 m en el multiespectral, con vista fuera del 
nadir de hasta 50 hacia cualquier dirección del azi- 
mul, ¡ta observar una misma área cada 2 0 
3 días. IKONOS orbita la tierra cada 98 minutos a 
una altitud aproximada de 400 millas. Con IKONOS 
se obtienen imágenes a 11 bíts (2048 niveles de gri- 
ses) lo cual genera un rango de información muy 
detallado. Space Imaging comercializa productos 
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con rango de de 4, 10, 12, 25 y 50 m 
A 
los efectos topográficos. 


PCIOrthoEngine 

Este software comercial usa el modelo desarro- 
llado por el Dr. Thierry Toutin del Canadian Center 
for Remote Sensing (CCRS), el cual fue probado 
usando otros dos métodos adicionales que son el 
método polinomial y racional polinomial, El prime- 
ro sólo cornge distorsiones planimétricas básicas 
y el segundo involucra la altura de los puntos de 
control para rectificar la imagen. Estos últimos no 
requieren información del sensor ni del satélite. En 
cambio, el método que involucra el modelo a util- 
zarse refleja la realidad física desde el punto de 
vista geométrico y corrige las deformaciones pro- 
ducidas por la plataforma, el sensor, la Tierra y 
algunas deformaciones producidas por la proyec» 
ción cartográfica, O sea, este último método consi- 
dera la información del satélite y del sensor para 
ontorectificación. 

La funcionalidad del software fué comprobada en 
una localidad de Toronto, Canadá con elevaciones 
comprendidas entre los 180 y 240 m. Allí fueron 
colectados 30 puntos de control distribuidos unifor- 
memente sobre la imagen obtenidos de ortofotos 
con 20 cm de precisión y un DEM a 2 m. Obtenien- 
do de la ortorectificación un RMS de 1,3 m y un 
error máximo de 3,0 m. (Toutin,2.000). 


3. LA METODOLOGÍA 


imagen IKONOS 

Las imágenes IKONOS (GEO, pan sharpened) uti- 
lizadas en este proyecto fueron adquiridas por el 
CPDI en el marco de un contratado por el 


bre de 1999. La escena cubre desde el sector Gra- 
moven al oeste, hasta el Country Club al este de la 
ciudad de Caracas, y desde el 23 de enero por el sur, 
hasta Galipán al Norte (véase Figura 1). La topogra- 
fía en el valle de Caracas varía desde los 890 m 
hasta los 1.200 m y en la zona montañosa sobrepa- 
sa los 2.000 msnm. El tipo de uso en la zona urbana 
comprende desde zonas residenciales unifamiliares 
de baja densidad hasta multifamikares de alta densi- 
dad altemando con zonas comerciales y administra- 
tivas. 

La imagen IKONOS fué adquirida el 18 de febre- 
ro de 2.000, El producto ha sido adquirido en for- 
mato GEO (pan sharpened) en proyección UTM 
19, datum WGS84 con una precisión nominal de 
50 m (CE90). Para la ortorectificación de la imagen 
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se usaron 50 puntos de control obtenidos por la téc- 
nica Diferencial GPS, un DEM con un tamaño del 
plxel de 20 m y el software PCI OrthoEngine. 

Se evaluó la precisión de la ortorectificación de 
la imagen IKONOS sobre un área urbana seleccio- 
nándose una subescena en la cual se muestra la 
distribución de los puntos de control (véase Figura 
5) que comprende desde las inmediaciones del 
Parque del Oeste en Catia hasta el Colegio La 
Salle al Este de la imagen IKONOS, para un total 
de 28 puntos de control y 5 puntos de chequeo 
sobre el área de estudio. 





Figura 1. imagen IKONOS A donde 
se representa hacia el sur, la zona Urbana de Cara- 
cas y al Norte el Parque Nacional El Ávila, Los 
GCPs están representados en Rojo y los de Che- 
queo en Cyan. 


Mediciones GPS 

El Sistema de Posicionamiento Global (GPS) 
está basado en la transmisión de información 
desde una consielación de satélites hasta los 
receptores GPS que están sobre la superficie 
lerrestre. La red del sistema GPS está dispuesta 
de tal forma que un receptor sobre la superficie 
terrestre puede captar un minimo de 4 satélites en 
cualquier momento del día y en cualquier parte de 
la superficie terrestre. Cada satélite describe órbi- 
tas circulares a 20.180 Km. de altura, con un perío- 
do aproximado de 12 horas. Como la altura de los 
satélites equivale a tres veces el radio de la tierra, 
cada receptor está en capacidad de recibir la señal, 
por un periodo aproximado de 6 a 7 horas. El uso 
de dos receptores grabando simultáneamente y 
captando la señal de casi todos los satéles es una 
manera efectiva de disminuir el error y aumentar la 
precisión del posicionamiento. 

En el presente trabajo fueron colectados 70 pun- 
tos GPS durante el mes de enero del 2.002. Las 
condiciones atmosféricas reinantes variaron entre 
nubladas y soleadas. Dos receptores GPS de uso 
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cartográfico fueron usados (GeoExplorer 3), Estos 
receptores garantizan precisiones de 20 cm cap- 
tando datos en modo Fase y corrigiendo diferen- 
cialmente. 

Previo a la Campaña GPS se definió como esta- 
ción base el instituto de Ingerveriía (FII), establecien- 
do sus Coordenadas Geodésicas a partir del punto 
IDEA de la Red GPS Metropolitana de Caracas de 
Orden C definida por el IGVSB. Mediante la técnica 
de recolección de datos en modo fase y el método 
de medición GPS Diferencial, se obtuvo asi el punto 
Fl con coordenadas N 1.151,274,30, E 731.444,13 
y 1.362,06 m sobre el elipsose con precisiones de 
20 cm en plansmetria y de 20 cm en altimetria referi- 
dos a la proyección UTM 19 datum La Canoa. Se 
venfico la precisión de los puntos de control reall- 
zando mediciones diferenciales a dos puntos del 
IGVSB: el IDEA a 1,24 Km y AVILA a 15,38 Km del 
IGVS8 obteniendo errores de 0,08 m en planimetría 
y de 0,03 m en altimetria, y de 0,06 m en planimetria 
y de 0,23 m en altimetria respectivamente ambos de 
la Red Metropolitana GPS de Caracas. 

Previo a las mediciones GPS se realizó la esco- 
gencia de los puntos de control (GCP') y de verífica- 
ción con alta calidad de reconocimiento en la ima- 
gen para que pudieran ser ubicados con la mejor 
precisión posible, (véase las Figuras 2, 3 y 4), evi- 
tando los lugares cubiertos por sombras producidas 
por edificaciones, así como las zonas no “vistas” por 
el sensor, las cuales impsden observar detalles de la 
imagen. Otro factor importante a considerar es evi- 
tar tomar puntos GPS donde se encuentren obs- 
trucciones verticales, ya que disminuye la captura 
de la señal desde los satélites. 





Figura 2: Punto de Control 304. Se puede Observar 
la posición precisa y nítida del punto sobre la ima- 
gen, que representa la esquina de una caminería de 
arcilla (so muestra dentro del circulo amarillo), 
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Para todos los puntos de control se tomó la 
siguiente información: Numero del punto, día, des- 
cripción del lugar, la hora de la captura, y la respeo- 
tiva minuta de campo. 

El procesamiento de los datos de la campaña 
GPS fué efectuado con el software GPS Pathfin- 
der Office 2.80 de Trimble. De los 70 puntos ubica- 
dos sobre la escena, 50 son puntos de control 
GCPs y 20 son puntos de verificación. Se presenta 
la precisión estimada por el software en la Tabla 1, 
la cual es de 1,22 m en planimetria y de 2,15 m en 
altimetría. 





Figura 3: Punto de Control 308, Plaza Candelaria. 
La marca oscura representa la estatua. La ubica- 
ción precisa del punto se observa en la esquina de 
una jardinera, notando la diferencia entre el jardin 
y el concreto, 





Figura 4. Punto 307. Colegio La Salle. Punto sobre 
al lindaro entre la cancha de Bólabol y la de Futbo- 
lito, y la Aw. Boyacá. 


Tabla 1: Precisión de los Puntos GPS. GCP puntos 
de control. GVP. Puntos de Verlficación. 
19408 los Ls A 
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Las alturas elipsoidales (h) fueron transforma- 
das a alturas ortoméincas (M) usando el Modelo 
Geopotencial Global EGM96, con una precisión 
nominal de 0,5 m, desarrollado por el National Ima- 
gina Mapping Agency (NIMA), el cual, realiza una 
corrección a la Ondulación Geoidal (N), cuando se 
le introduce la altura elipsoidal al programa. (NIMA 
EGM96 Calculator Ver 1.1). 


DEM (Digital Elevation Model) 

El DEM utilizado fué producido a partir de las cur- 
vas de nivel con intervalos de 40 m de tas cartas a 
escala 1:100.000 editadas por el IGVSB. El mode- 
lo fué construido con un espaciamiento en la cua- 
dricula de cada 20 m para eliminar los errores de 
interpolación en ta imagen con las curvas de nivel. 
Las características del DEM son las siguientes: 

La precisión planimétrica en el 90% del área es 
de 50 m, 

La precisión altimétrica es de 20 m en el 90% de 
los Casos. 

El operador digitaliza con un error de 0.1 mm lle- 
gando inciertamente a 10 m. para escalas 
1:100.000. 

En el proceso de interpolación para generar el 
DEM resultó dificil cuantificar la contribución del 
error, la cual fué estimada como irrelevante. 

La topografía que caracteriza al sector del Valle 
de Caracas que comprende el área de estudio, 
está compuesta por pendientes medias a altas 
hacia el norte y noroeste, donde se localizan las fal 
das de la Serranía del Ávila. Esta tiende a ser más 
plana hacia el este, a medida que aumenta la dis- 
tancia del pie de monte. Hacia el sur y suroeste, se 
localizan las colinas donde se emplaza la urbani- 
zación 23 de Enero, lo cual se traduce en pendien- 
tes medias a altas. Por otra parte, dentro del paisa- 
je resaltan algunos quiebres de pendiente de gran 
relevancia como el Cañón de Tacagua y los cau- 
ces de taz Ovebradaz Manicomio, Lídica y Catu- 
che. 


4, RESULTADOS 


La imagen orlocarregida se puede observar en 
a Figura 5, yse sesume enla Tahla 2 Donde enn la 
arto- ectificación se obtuwa un errar RAMIS an de 
4,50 pixales y en Y de 2 53 pixales y al de los pun- 
los de chequeo en X as de 27,69 pixeles y er Y de 
12.06 pixe es. 

Para mejorar la distribuciór espacial del error se 
seleccionaron los puntos de control con mayar pre- 
sisión. Se escogió como tácnica de-replanteo del 
ohrel el método bilineal, porque no se observaron 
distorsiones visibles de líneas a alta resolución. Se 
utikzó como datum la Canca. 
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Para la verificación de la imagen IKONOS orto- 
rectificada, se utilizó: puntos de chequeo GPS, las 
imágenes ortocorregidas en el CPDI (Ruiz, 2.000) 
y los ortofomapas (IGVSB), 

Los Onrtofotomapas del proyecto Cartocentro 
del IGVSB, levantado a un pixel de 2,5 m con una 
precisión planimétrica de 10 m (excelente correta- 
ción visual con los puntos GPSs en toda el área) 
con los que fueron usados como referencia para 
observar la relación existente con la imagen orto- 
corregida, donde había una coincidencia en la 
mayor parte de la zona de estudio (ver Figura 6), 
exceptuando al Este del sector Ruperto Lugo que 
se localiza en la zona Norte de Los Fralles (ver 
Figura 7) y una pequeña variación en la parte Sur 
de la imagen especificamente en la zona del 23 
de Enero. 





Figura 5: imagen Ortocorregida representada 
sobre un Modelo Digital del Terreno, 





Figura 6: Comparación entre la imagen Ortocorre- 
gida (Vista derecha) y los Ortofotomapas (vista 
izquierda). 
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Figura 7: Comparación entre la imagen Ortocorre- 
gida (Vista superior) y los Ortofotomapas (vista 
inferior). 


Comparación con la imagen ortocorregida por 
el CPDI (2.000), con la imagen ortocorregida en 
este trabajo. (Ver Figura 8). En general se puede 
observar una buena correlación entre las dos imá- 
genes ortocorregidas exceptuando en la zona 
Norte de Los Frailes como se mencionó anterior- 
mente, 





Figura 8: Comparación entre la imagen 
Ortocorregida en esto trabajo (vista derecha) 
y Ortocorregida por el CPDÍ en el año 2.000 
(vista izquierda). 


Los puntos de Chequeo GPS fueron converti- 
dos a una capa vector para luego ser correlacio- 
nados conlas posiciones que deberian llevar ori- 
ginalmente, cuando fueron medidos. Mostrán- 
dose errores totalmente distintos alos dados por 
el RMS calculado por OrthoEngine. (véase la 
Tabla 2). 
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Tabla 2: Precisión de la Ortorectificación. Dif es la 
diferencia sobre la imagen corregida de los 
puntos de Chequeo. Prec DGPS es la precisión de 
la corrección diferencial. RMS Ortho es el error 
RMS de la Orto rectificación. 
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Análisis de Resultados 

Para la ortocorrección de la imagen del área 
urbana de Caracas se presentó varias condiciones 
negativas, en cuanto a que la imagen no posee 
parámetros orbitales, tiene un ángulo de 
incidencia de la imagen de 50" y un DEM de baja 
resolución para esta zona por poseer una 
topografía vanada, 

Sin embargo, cuando se realizó la 
comparación visual de la imagen IKONOS con 
los Ortofotomapas del proyecto Cartocentro, en 
la cual se observó una buena similitud entre los 
elementos espaciales sobre casi toda la imagen, 
exceptuando la zona norte el sector llamado Los 
Frailes donde existen desplazamientos en 
dirección Norte y Este del orden de 15 m lineales 
aproximadamente. También existen 
desplazamientos hacia el sur de la imagen en la 
zona del Barrio Santa Clara, de 
aproximadamente 17 m., Si observamos la 
Figura 9, el DEM de la zona permite definir que 
hay cambios de pendiente en la topografía lo que 
definitivamente afecta la ortocorrección de la 
imagen. 

Cuando se comparan los ortofotomapas y la 
imagen ortocorregida con la cartografía digital, se 
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observaron diferencias infenores a los 2 metros 
debido a que en esta zona la topografía no presenta 
cambios bruscos de pendiente. Como sucede en 
las inmediaciones del Parque del Oeste en Catia, 
en el Barro La Linea, Los Flores, Gato Negro, 
Lídice, Manicomio, algunas zonas del 23 de Enero, 
La Pastora, Altagracia, San José, La Candelaria, 
Maripérez, San Bernardino, Sarría, Guaicaipuro, A 
medida que se acerca al pie de monte del Parque 
Nacional el Ávila tas diferencias se incrementan 
sobre todo en el área de Puerta de Caracas donde 
se identifican pequeñas deformaciones en 
dirección Este cuando se observan dos laderas, 
este fenómeno es muy visible en los barrios 
Macayapa y Resplandor, encontrándose estos 
sobre el parque nacional el Ávila 





Figura 9: DEM de la zona de estudio. Se observan 
los cambios de pendiente existentes en las zonas 
problema, que se presentan en la parte inferior y 
superior izquierda de la imagen. 


En cuanto a los puntos de chequeo el 005 y 013, 
el primero ubicado en la calle principal de Mirador 
y el segundo en los altos de Lídice. Estos poseen 
diferencias que están en el orden de los 2 m, el 
109 ubicado en La Candelaria (centro de 
Caracas) tiene una diferencia de 2,16 m, 
Adicionalmente se midieron otros tres puntos 
ubicados a una distancia aproximada de 100 m, 
donde se hace mención al punto del Hotel Ávila 
que tiene una diferencia por debajo de los 2 m y 
los otros dos por estar muy cercano al Ávila tiene 
una precisión de 2,6 y 2,7 m, con esto se 
evidencia lo importante de un DEM preciso 
cuando se tiene una superficie topográfica muy 
variada, tomando como caso especial al colegio 
La Salle, que al estar muy cerca al Parque 
Nacional el Ávila y ser el punto de control extremo 
a la escena, tiende a producir deformaciones por 
la falta de un DEM preciso. 
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5. CONCLUSIONES Y RECOMENDACIONES 


Para la imagen IKONOS (GEO, pan sharpened) 
realizada en el área urbana de Caracas, NO se 
obtuvo la exactitud cartográfica requerida para la 
escala 1:5,000, porque el RMS de los puntos de 
control fué de 4,90 m. A pesar de eso, al comparar 
las zonas planas con los Puntos de Verificación y 
con otros productos ortocorregidos, se concluye 
que ésta alcanza las precisiones de un producto a 
escala 1:5.000 para zonas con cambios de pen- 
dientes inferiores al 8%. 

Se recomienda el uso de Imágenes Orthokit 


para trabajar enáreas urbanas. 
Se recomienda el uso de DEM's de mayor reso- 
lución para ortocorregir ¡ Ikonos donde 


existan cambios evidentes de pendiente. También 
$e recomienda evaluar la posibilidad de extraerlos 
a partir de estéreo imágenes de sensores como el 
LIDAR, SPOT e KKONOS. 
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